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I. INTRODUCTION 

Modern control methods allow the control engineer to 
develop a control system to regulate every parameter ofa 
physical systen. He has been aided by computer simulations 
of non-linear systems and linear approximations of those 
systems. The desired contrcl system can be developed 
parallel to or even ahead of the physical system which is to 
be controlled. 

In .many design applications areas, the engineer 1S 
confronted with designing a ccntrol system which can be 
represented analytically by a very large number of first 
order, linear, time invariant differential equations. This 
is especially true with highly complex and aero-elastic 
alrcraft and space vehicles. The analytical model is a 
prime ingredient in the controller design process for any 
design technigue. 

Many practical limitations exist, such as computational 
requirements, that require the order (number of equations) 
of this complex model to be reduced for analysis, synthesis 
and actual implementation of the control system. Some method 
is reguired to analyze these dynamical models and establish 
Simpler systems which include elements critical to the 
desired control function. Weehout such @=yslmpsal i leat. on, 
application of design procedures can result in highly 
complex and parameter-sensitive controlled systems. With a 
reduced order model, regulator synthesis procedures become 
more intuitive and far less sensitive to parameter varia- 
TOT) « The reduced order model is far easier to handle in 
formulating the control system but must characterize the 
physical system with sufficient accuracy such that perforn- 
ance objectives for the controlled physical system can be 


met by designing control laws from the reduced order model. 
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The intent of this thesis is to present the mathematical 
basis for creating such a reduced order model and to develop 
the actual computer software which creates that reduced 
model based on a given large order systen. This computer 
program (OPTRED) is developed to interface with an existing 
control analysis FORTRAN prograf named OPTSYS. This thesis 
will discuss the modifications or additions to the OPTSYS 
program which are necessary to that interface. 

The OPTimal SYStems control program was originally 
developed by Hall [Ref. 1] for the study and application of 
optimal systems control theory. Later modifications tc the 
program were made by Walker [Ref. 2] and Liu [Ref. 3] of 
Stanford University. Hoden [ Ref. 4] modified the OPTSYS 
program to present a user-friendly, interactive version of 
the program (OPTSYSX). Diel [Ref. 5] introduced changes 
which enabled the user to save various data for re-entry 
into the program and also to create data files which were 
necessary for the execution of atime response program 
(OPTCALC). Further modificaticns have been made by Laptas 
[Ref. 6] which were necessary tc create data sets for input 
to his OPTGRAPH progran. This FORTRAN program enabled the 
user to obtain Pole-Zero, Root-Locus, Nyquist, Bode and 
Nichols flots. 

It is assumed that the user has a basic understanding of 
the fundamental concepts of control theory and . optimal 
systems design. The symbolcgy conventions of Bryson 
[Re£f. 7] are used in the discussion of program operation and 
system descriptions. 

An overview of the OPTSYSX prograna, 1ts capabilities, 
and the modifications to that program which are necessary 
for the interface with OPTRED is presented first. This is 
followed by the mathematical kasis upon which OPTRED is 
formulated and by a full description of the operation of 


that. program. Finally, examples of large order control 
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systems and the corresponding reduced modeis which OPTRED 
generates are presented and analyzed. Those systems include 
a sixteenth order model of the F100 turbofan engine anda 
Ninety-eighth order model of the X-29A aircraft. Complete 
program listings for the OPTSYSXK program and the OPTRED 


program are included as Appendices A and B respectively. 
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IIT. THE OPTSYSX CCUPUTER Pica 


A. PROGRAM OVERVIEW 


OFTSYSX is a double precision, interactive FORTRAN 
program which employs modern ccntrol theory analysis tech- 
nigues. It is developed to be compiled and executed by the 
Naval Postgraduate School's IBM 3033 System 360/370. Its 
primary capabilities include the calculation of the open 
loop eigensysten, and the fixed closed loop system; the 
synthesis of optimal reguiators or filters; the power 
spectral density, and modal distribution computations. 

The fundamental system equations used by the OPTSYSX 
program for its computations are of the state variable forn. 


The basic system equations are: 


system model 

x = [F]*x + [G]*u + [ GAM ]*w (2245 
output equation 

Y o> GPa ee (2-2 
measurement eguaticn 

z= [H]*x + v (2.3) 
estimator equation 


xe = [F]*xe + [G]*u + [K]*(z - [H]*xe) (2.4) 
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open loop transfer function 

Case P=? *, S (2.5) 
closed loop noise transfer function 

figs (Siilg-— (eget ican | (2.6) 


compensator transfer function from measurement to input 


PemiSset PLE LG te 14K eth )[=***] (2.7) 
where 
x = state vector (Ns X 1) 
x = derivative of the state vector (Ns X 1) 
u =—GOMmenoL Vector (Ne X 1) 
y = output vector (No X 1) 
Zz = measurement vector (No X 1) 
W = white process noise vector (Ng X 1) 
[ F ] = open-loop dynamics or plant matrix (Ns X Ns) 
[G] = control distributicn matrix (Ns X Nc) 
{Gam} = state disturbance distribution matrix (Ns X Ng) 
ae = measurement distribution matrix (No X Ns) 
[D] = control feed-forward matrix (No X Nc) 
irc } = control gain matrix (Nc X Ns) 
et ] = identity matrix (Ns X Ns) 


OPTSYSX 1s an extremely large and complex program which 
contains over 4000 lines of code. Its existing standard 
dimensions will accomodate a thirty-second order system with 
up to 10 controls and 10 observations or measurements. LO 
execute even this version of OFTSYSX the user requires one 


megabyte of virtual storage. Any further significant 
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additions to that program's size would cause user difficulty 
in obtaining sufficient memory for program execution. It is 
for this reason that Dielts OPTICALC, Laptas' OPTGRAPH and 
the program developed by this thesis have been developed as 
independent programs which depend upon various data created 
by OPTSisk. The minor modifications this thesis introduces 
to the OPTSYSX program are not detrimental to its original 


capabilities. 


Be OFPTSYSX MODIFICATIONS 


This thesis is ccenfronted with the problem of analyzing 
very large control system nodels. Therefore, the additions 
of Hoden'ts subroutine SETUP and Diel's subroutines RDMATF, 
RDMAT, and WRIMAT were crucial to system analysis. These 
routines make multiple computer analysis runs’ possible 
without painstaking re-entry of system parameters’ and 
matrices. In particular, the system "F", "“G", ij" eae 
"Gam" are some or the matrices which are stored for re-entry 
into OPTSYSX or stored as data for use by other analysis 
programs. Subroutine WRIMAT is responsible for the storage 
of these matrices inthe data file "OPTMAT". It is this 
file which OPTRED utilizes tc input all pertinent full 
system data. 

In many control applications, a control feed-forward 
matrix ("D") as present in the system model. Therefore, 
this thesis makes proviSions within OPTSYSX which enable the 
user to save the "D" matrix in the same manner by which the 
other system matrices are saved. Subroutine SETUP, which 
enables the user to input system matrices from file data, 
was Changed to enable the user to input the "D" matrix. 
This consisted of adding "D" to the routine's parameter list 
and also added a read statement to input that matrix frona 


data file. Present routines in OPYSYSX which "save" the 


abe 


== 


uP ong 6mm, 6©6Uacnd )6©"GAM" matrices were complemented by the 
addition of routines which also save the "D" matrix. These 
new routines are analagous to existing "Save" routines. In 
addition to these routines a new flag, ISAD, was added to 
eye Existing flags ISAF, ISAG, ISAH, and IGAM. This new 
flag provides the same function in "Saving" matrices as do 
the existing flags. For program continuity, it was also 
necessary to add ISAD to subrcutine INNER, which is the 
driving subroutine in the OPTSYSX progran. 

Subroutines RDMATF, RDMAT, and WRIMAT also required 
minor changes to accomodate the "D" matrix. These changes 
consisted of the addition of "D" to the subroutine parameter 
lists andthe addition of read or write statements for 
fe~ner re-entry of "D" into OPISYSX or for writing "D" to 
the OPTMAT data file. The existing flag, IFDFW, was also 
added to the parameter lists of these routines. This flag 
indicates whether or not a "D" matrix will be input by the 
user. OPTRED uses this flag in its computational methods in 
the same manner as does OPTSYSX. Appendix A lists’ the 


OPTSYSX program as modified by this thesis. 
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III. ANALYTICAL DEVELOPMENT OF A REDUCED ORDER MODEL 

There are numerous methods of order reduction for 
control systems as discussed by Enns [{Ref. 8] of Stanford 
University. This thesis uses the theory discussed by DéeHoff 
and Hall [Ref. 9] which presents a method of reducing the 
order of the linear system based on dominant states chosen 
by the designer. This dominance is determined by modal 
analysis of the eigensystem under investigation. Dougherty 
[Ref. 10] utilized this same methodology in the development 
of a computational model used in the analysis of the F100 
turbofan engine. A summary of this technigue 1s provided 
below. 

The method of reduction involves finding the eigenvector 


transform matrix, T, such that 
eee (3 49) 


where A is the diagonal eigenvalue matrix. Then, by 


defining an alternate state vector, Z, where 
X = TZ (3225 


the original system may be re-written in modal coordinates 


as 


Z = A¥Z + T-1*G*U (3.3) 
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The eigenvalues which are tc be included in the initial 
reduced order model are determined through modal analysis. 
The eigenvalue matrix can be re-ordered by constructing a 
matrix Q, which has the value 1 in the position corre- 
sponding to the eigenvalue to ke included (column) and the 
re-ordered position (row) of that eigenvalue. FOr example, 
if eigenvalue 20 of the original A matrix was to be 
re-ordered to the first position, a 1 would be placed in the 


(1,20) position of the Q matrix. Now, the reordered matrix 


At = ¢C*A (3.4) 


is incorporated in the linear model 


O#Z = A'#Z + Q¥T-14G*U (3.5) 
letting 
Z' = QxZ (3.6) 
thus 
ZY = AY*kQ-14Z" + Q*T-14GKU 1259) 
and 
X = TeQ-1«z! (3.8) 


A matrix R can be constructed such that 
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1 
X* = R*K = TZ! = |--=- (3 oe 
X 
Za 
This X* vector is the reordered state vector, where the 


elements contained in the subvector x, are those states 
associated with the designers choice of eigenvalues while 


the subvector X.~. contains all other states. 


2 
The linear system may now be rewritten 


~ 


X' = R*¥F*R-1¥*X' + R¥G¥*O (3.10) 
and 
Y = H¥R-1*X' + D¥*J (3.15 
Now define 
Ft = R*¥F¥R-1 (3.72 
G' = R*G (3.73 
H' = H*¥R-1 (3.14) 
and 
De 
Pl eee 
Fr = amen aa (3. 13 
. 2 


G*t = —— = 


Pieene states in xX, truly )model the full order 


then the subvector Xx. will be essentially zero and 


= x + x + * 
= Pu a P12 * 5 a : 


= * + a + * 
: ol a} P22 a 25 ° 


Solving for X5 in Equation 3.19, one has 


= - -1x * —-1*G x 
P.5 Poi Xx + P55 Dy U 


This 1s incorporated into Equation 3.18 and 


ee (hor ELE Monee Corm= ee oe te ) Lf 
1! ( eel 29) 22 21) al ( 2 a) 


and the output equation becomes 


Y = (H.-H _*F -1%F + M( DoH oe FE 1k 
Ee OD et ee pn 


The reduced order model can now be written as 


X. = Fr*¥xk + Sr*¥vJ 
IL A 


Za 


(3.16) 


(3.17) 


SyYSten, 


SiG 118) 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


Y = Hr*xX 
alt 


where 


Pr--= 


HE 


I) 
a8 
i 


+ Dr*U (3.24) 
yes 2 (Seay 
* —1* 
12 2m ane (332 
* -~1x 
Ho*F,.-t*F,, (3.29) 
-1* : 
5-1 *6, (3.28) 


= - H * 
Dr D PF, 


This is the mathematical prccedure which OPTRED uses in 


computing reduced order models. 


is included as Appendix B. 
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The complete OPTRED listing 


A. PROGRAM OVERVIEW 


OPTRED iS an interactive, double precision FORTRAN 
program which uses a specifically formatted data file to 
input full system data. The required matrix inputs from 
this file are the full systen "F", Git and. ps 
Matrices. The system flags and parameters which OPTRED must 
also read for program execution are Ns, Nc, No, and IFDFYW. 
The OPTSYSX program creates OPTMAT data which is the source 
of this required data. Through interactive input, the user 
specifies all other data which is required to create reduced 
model data. 

OPTRED also reads from the CPTMAT file the noise, § feeéed- 
back control and the output and control cost weighting 
Matrices. These are not reguired inputs but only input to 
the program to facilitate creation of a reduced order data 
file which is identical in format to the OPTMAT data file. 
The final result of OPTRED is the creation of OPTMATR, a 
data file which contains the reduced system "F", "G", "HN" 
and "D" matrices andis formatted for re-entry into the 
SEioYSX progran. 

It iS possible for the user to construct a data file 
hamed OPTMAT to provide input data to OPTRED. However, only 
in rare applications would the full system not be analyzed 
prior to the development of the reduced nodel. OPTSYSX is 
therefore the simplest and most preferred method of data 
generation. 

to0.)6r execute the OPTRED prcgran, a Minimum of three 


separate files must be immediately accessible to the user: 
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1) The compiled text file of OPTRED 
2) The OPTMAT data file 
3) The OPTRED EXEC file 
The first two required files are self explanatory. The 
OPTRED EXEC file sets up the required read/write file defi- 
nitions and calls the necessary library functions for execu- 
tion using the IBM 3033 Systen. This EXEC allows the user 
to execute the OPTRED program by Simply entering the werd 
NOPTRED" at thes cerm incr The program then continues in a 
user-interactive mode until completion as described in the 
next section. The OPTRED EXEC file is listed in Appendix C. 


B. PROGRAM OPERATION 


OPTRED waS written to satisfy two specific objectives. 
The first is to provide a methcd to assist in the analysis 
of actual large order control systems. The second is to 
create a program to be used aS an instructional tool for 
Students studying controls related problens. 

To accomplish these objectives the emphasis while 
creating OPTRED was to make the program as user oriented as 
possible by eliminating data entry ambiguities and providing 
features which prevent the user from other inadvertant or 
invalid program entries. A serious effort was also made to 
Minimize the amount of virtual memory required to execute 


thUS ~Eeogran. 


1. Program Composition and Flow 


OPTRED 1S composed of one Main program and = seven 
subroutines. The program can be divided into three fasic 
catagories: 

1) File Data Input 
2) Interactive Data Input 
3) Cal Giiageron 
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Both file and interactive data entry is accomplished 
by the main program and by the primary subroutine, REDUCX. 
The actual computation of the reduced order model is 
performed by subroutine REDUCX which utilizes subroutines MI 
and MAMULT for its matrix inversion and multiplication 
requirements. The remaining Subroutines are utilized 
throughout OPTRED for interactive entry of numerical -and 
character string data. The f~rogram requires no external 
subroutines for its execution. 

The main program begins by presenting the user with 
a brief operational description of OPTRED. himenas poLnt, 
the user is given the option of having the program's general 
results repeated to the screen cr printed to a iisting file. 
OPTRED creates a permanent data file which contains the 
reduced data regardless of the user's option here. The 
program then reads full system flags and parameters, 
excluding matrices, from the OCPTMAT data file. A brief 
description of the full system under investigation is then 
presented as an indication to the user that the desired data 
is being input. This description includes the order of the 
system, the number of controls and observations and whether 
or not a systen "D" matrix will be input. The latter is 
indicated either by " "or "nct" in the following example. 
The program then prempts the user for entry of the desired 
order of the reduced systen. 

The main prcgram presentation as it appears during 
execution: 

OPTRED 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SYSTEM 
og GY HY, AND "D™ (IF INPOT) MATRICES. 


THE FULL SYSTEM MATRICES MUST BE READ FROM A FILE NAMED 
"“OPTMAT DATA" AS CREATED BY THE OPTSYSX PROGRAM. 


YOU MUST ALSO ENTER THE DESIRED REDUCED ORDER (NUMBER OF STATES) 
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AND THE ACTUAL STATE #*"S (IN ASCENDING ORDER) WHICH REPRESENT 
THE REDUCED MODEL. THE ORDER OF THE REDUCED ORDER MODEL MUST 
BE LESS THAN THE ORDER OF THE FULL SYSTEM. 


DO YOU WISH TO CONTINUE? 
TYPE “YES™OOKS Haye 
DO YOU WISH RESULTS TO SCREEN O8 DISK? 


HOTE: A DATA FILE CONTAINING THE REDUCED ORDER DATA 
WILL BE GENERATED REGARILESS OF YOUR ANSWER. 


TYPE "S" FOR SCREEN OR "OD POhe ce 


THE ORDER OF THE FULL SiStiawts. 
THE NUMBER OF CONTROLS IS5:3 

THE NUMBER OF OBSERVATIONS IS: 

A "D" MATRIX WILL BE INPOT. 


DO YOU WISH TO CONTINUE? 
TYPE CUYES™ 7OR> NON. 


ENTER THE DESIRED REDUCED ORDER OF THE "F" MATRIX. 


Subroutine REDUCX is now called and immediately 
reads the full system matrices from file data. The user is 


then prompted for entry of the significant states which will 


compose his reduced model. This program makes a provision 
for entry of Significant states  Sfn0n a specifica 
formatted, user created data file. At this point ine 
program heis given that option. Once the significant 


states have been entered, the program repeats those states 
to the screen and the user is given the option to make 
changes at that time. If the significant states are satis- 
factory to the user, the program continues with the computa- 
tion of the reduced order model. At the completion of 


OPTRED, the reduced system data is created and can fe exan- 
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Pe CecOmmstantl ity, sconutrollalility and observability by 
other analysis technigues as discussed in the following 
section. Prior to exiting the OPTRED progran, the user is 
Meevaded the Option Of tre-executing OPTRED if he finds 
initial results unsatisfactory. A typical recording session 
mireh depicts this program ficw is presented in the next 


section. 
2. User Protection Features 


The user of OPTRED is protected from abnornal 
program termination in several ways. Initially, the progran 
presents its capabilities and user required entries. If he 
is not prepared to continue at this point, the user is 
offered the option of exiting the progran. Next, the full 
system parameters are presented and the user is once again 
able to exit if these parameters are unexpected. 

OPTRED reguires that the order of the reduced systen 


ke greater than zerc but less than the order of the full 


systen. If the user inadvertantly enters a number which is 
out of this range, the program issues a warning to that 
effect and allows recovery. The computation of the reduced 


order model relies upon the construction of a reduced state 
matrix. OPTRED requires that these states be entered in the 
order that they would normally appear in this matrix. For 
this reason, the user must enter the significant state 
humbers in ascending crder and those state numbers must ilie 
in the range from one to the rumber of states in the full 
systen. If the user inadvertantly enters these _ state 
humbers out of order or out of range the program frompts a 
warning and allows recovery. 

The method of reduced order computation involves the 


Qnversion of the state submatrix as described in 


Foe 
22 
mrapter TiT, Equations 3.25 threugh 3.28. The nature of the 


plant matrix and the selection of reduced order states can 
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be coupled such that this matrix wall be singular ODwieae 
invertable. OPTRED will detect this condition and issue 
information to the user that a reduced order model cannot be 
calculated for that full system from the given desired 
reduced states. Reduced data will still be computed but the 
user is warned that this data is invalid. 

Subroutines RDINT, RDREAL and RDCHAR are responsible 
for the interactive input of data and expect an integer 
input, a real number input and a logical "YES" on "NO" 
respectively. If the user inadvertantly makes an incompat- 
ible entry these subroutines issue warnings and allow 
another opportunity for data entry. The entry of a null 
line is included in these improper actions and the entry of 
two consecutive null strings will cause termination of the 
program. This function allows the user a further means of 


exiting OPTRED if he so desires. 
3. Large Order Systems 


The analysis of large crder systems, particularly 
the ninety-eigth order model, presents a major problem to 
the user with limited virtual memory assets. The longitu- 
dinal control system analyzed by this thesis has a (98 X 98) 
"NFW matrix, a (2 X 98) "H" matrix and a (98 X 1) "G" tate 
To analyze this systen, OPTSYSX must be executed in its 
"increased dimension" form. To compile and execute this 
version of OPTSYSX requires 2.5 Megabytes of virtual memory. 

Both OPTRED and OPTGRAPH are also dimensioned to 
enable analysis of a ninety-eigth order model. When 
OPTSYSX, OPTRED, OPTGRAPH and their peripheral data are all 
used during a typical analysis session, very careful atten- 
tion to file and virtual memory management 1s mandatory to 
prevent the user from exceeding his virtual machine's 


Capacaey. 
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C. INTERFACE WITH OTHER ANALYSIS PROGRAMS 


To sufficiently analyze a ccntrol system it is desirable 
me, utilize OPTSYSX, OPTCALC, CPTGRAPH and OPTRED. It may 
also be desirable to conduct this analysis during a Single, 
continuous interactive computer Session. However, in the 
analysis of very large systems this concept presents serious 
difficulties in the areas of program flow, data management 
and physical memory assets. For this reason, OPTRED has not 
been developed to be fully autcmatic in its interface with 
these other analysis programs. 

The OPTRED program requires the existance of the OPTMAT 
data file. AS previously discussed, OPTRED generates the 
OPTMATR data file which contains the reduced model data. To 
analyze this data using OPTSYSX, the user must now rename 
the OFTMATR data file for entry into OPTSYSX as OPTIMAT data. 
Care must be exercised at this point to either erase the 
full system data file or to rename it to prevent ambiguous 
data files. Although this procedure may seen cumbersome at 
mIEst, it 1S easily accomplished in the XEDIT mode and 
provides the user with the capability of maintaining several 
System data files which are readily available for analysis. 

After OPTSYSX processes the reduced data, additional 
data files are created which enable the execution of other 
system analysis programs. OPTGRAPH and OPTCALC utilize the 
OPTNAT and OPTGROL data files, respectively, to provide the 
functions described in chapter 2. When the user is investi- 
gating several systems it is both prudent and convenient to 
rename these files after their use for the purpose of future 


analysis. 


De. EXAMPLES OF ORDER REDUCTION 


Two large order systems will be presented to evaluate 


reduced models generated by OPTRED. The first system iS a 


Zz 9 


Sixteenth order model of the F100 turbofan engine. The 
final example 1S a hninety- eighth order model of the X-29A 
aircraft's longitudinal control systen. Bode analysis was 
conducted for the full and reduced models and graphical 
comparisons follow at the end of the chapter. Recorded 
terminal sessions for each example will be presented to 


fully illustrate actual program operarcrom. 


ey) 


1. Example of a Sixteenth Order System 





The F100 turbofan engine was chosen for initial 
analysis and [Ref. 9] describes this system in detail. The 
system "F'>. (Gl, "HY and "D" matrices were obtained from 
rRef. 9: pp. 83-85] and the method of data entry to OP tea 
is depicted in subroutine SETUF found in Appendix A. The 
selection of Significant states is based upon a desired 
control bandwidth of 1 to 10 Hertz and following modal anal- 
ysis these states were chosen as 1, 2, 5, 11 and 16. The 
following 1S a computer terminal system in which a fifth 


order model is generated. 


BEGIN RECORDING OF TERMINAL SESSION 
Re T=050170-03 220-0 5-4c¢ 
OPTRED 
FILEDEF 05 TERM 
FILEDEF 03 DISK REDUCI DATA A1 
FILEDEF 06 DISK OPTRED LISTING A1 
FILEDEF 07 DISK STATES DATA A1 
FILEDEF 09 DISK OPTMAT DATA A1 
FILEDEF 10 DISK OPTMATR DATA A1 
GLOBAL TXTLIB FORTMOD2 MOD2EEH IMSLDP NONIMSL 
LOAD OPTRED ( START 
EXECUTION BEGINS... 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SY oa 
NE WGN MA AND "DU (IF TN 2Ul ee eee 
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aul ot olteM MATRICES MUST BE READ FROM A FILE NAMED 
POPTMAT DATA™ AS CREATED BY THE OPTSYSX PROGRAA. 


mee MUS? ALSO ENTER THE DESIRED REDUCED ORDER (NUMBER OF STATES) 
mip THE ACTUAL STATE #'S (ZN ASCENDING ORDER) WHICH REPRESENT 
me neELUCED MODEL. IHE ORDER OF THE REDUCED MODEL MUST BE 

meoo THAN THE ORDER OF THE FULL SYSTEM. 


WOmrOU Wo TO CONTINUE? 


ioe fh” OR, NON. 


YES 
VOmrOUMmNionl REoULIo ITO SCREEN OK DISK? 


NOTE: A DATA FILE CONTAINING THE REDUCED ORDER DATA 
WILL BE GENERATED REGARDLESS OF YOUR ANSWER. 


ite wrOneoCh Eh EN CR "D" FOR DISK. 


$ 
Pmonvnm OF THE PULL SYSTEM IS: 16 
HMieNGMeER Or CONTROLS IS: 5 
POE SNUB ER OFS OBSERVATIONS ISs 7 
he DM MATRIX Wil. jo Ome 6. B20 aie 
OO YOU WISH TO CONTINUE? 
YE Dey Bow. OR ior, 
nH ES 
Tier eliotn REDUCEDSORDER OF THE "FEF" MATRIX. 
es 
18 


*¥e** WARNING: REDUCED ORDER MUST BE GREATER THAN 0 ****x 
AND LESS THAN 16 


ete litter ork De REDUCED ORDER OF THE "F™ MATRIX. 
es 
e) 
meee. OU WISH TO INPUT DESIRED STATES FOR YOUR REDUCED ORDER 
MODEL FROM A DATA FILE? 
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DATA FILE MUST BE NAMED "STATES DATA Al IN FIXED FORRae 
THE READ FORMAT IS "1315" PER 72 CHARACTER 2 ae 
TYPE "YES" OR “NOE 
NO 


ENTER THE "N" SIGNIFICANT STATES WHICH REPRESENT THE Regge 
MODEL. ENTER STATE #°S IN VASCENTING Clee 


STATE # 1 = 
i 
1 
STATE # 2 = 
63 
Z 
STATE # 3 = 
2 
3 
STATE # 4 = 
2. 
11 
STATE # 5 = 
G 
16 
THE REDUCED MODEL STALIES “Sie: 
1 Z 3 ile 16 
DO YOU WISH TO CHANGE ANY OF THE SIGNIFICANT Sie 
EYP GOV Y ES OR ome 
YES 
ENTER THE N-TH STATE @ 20) BERSHIGEoe 
ts 
3 
ENTER NEW STATE # 3 
5 
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Pion eouceu MODEL STATES ARES 


1 Z 


5 11 


16 


DemrOouentoH TO CHANGE ANY OF THE SIGNIFICANT STATES? 


NO 


iYP Er 


ny ES ij OR uno". 


eouUTPULT 10 THE SCREEN FOLLOWS) 


MiteveolreD REDUCEDZORDER 2S: 95 


Ths REDUCED NODEL STATES ARE: 


1 Z 


iit REDUCED 


aoe V209D+00 
fez 243D-01 
mee? 200D-03 
4.1314D+00 
me 0? JOD+00 


fis RELUCED 


—-3.5445D-02 
—4.4094D-01 
fee > 26D-02 
fe SO2D+01 
pale C909D+00 


THE REDUCED 


io 231D-01 

6.-4607D-03 

0 

m2 35D-05 
ereoo /0D-05 

feo 598D-05 

fa2415D-06 


5 11 


2-4726D+00 
-1.6271D+00 
22055 1De03 
-5.3854D+00 
262520 D+00 


-1.4080D+02 
aie hee D+ 0) | 
-3.7646D+01 
Seale ot OZ 
a2 OO DOG 


1.0750D-0 1 
—9e9442D-07 

0.0 

i551 2D>05 

1.4796 D-04 
ao. 0S D= OO 
Osmo D—O7 


16 


Pipa Ne MAC RX re?) 


to 


=e oot OZ 
4.1060D+t01 
-7.4276D+00 
aco) Dt 0-5 
1.14€7D+03 


CONTROL DISTRIBUTION MATRIX 


res Ga) | 
5.5 06D 00 
V2 ob 04 
4.9385D+01 
-3.0436D+01 


lo 7 Some 0 2 
Seow p— 0S 


0.0 
SAG ee sD 6 


ha TeOn 03 
mete! DOs 2 


ee Or oo DU 2 
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S10 9 13 D+00 
4.2898D+00 
=a O95 D= 02 
aS OO 
1.9087D+01 


("G") IS: 


Zea oT DEO 
Sale OOO DO 
= Te oweo) — 02 
= 7.00 15D 0 | 

Eyooor pt 0 | 


CUP U Oot RE BUTTON MALT RIX (“H™) IS; 


1.2088D+00 
SA (0) 1 ZA NS 

1.0000D+00 
-1.3900D-904 
oso O05 
es 0G 0—O4 
-4. 9442D-05 


Shs USS A Oe 
a2 022o00—0 1 
=9..45Z24D—03 
-2.7911D+00 
= Woo 71 D+0] 


Saison opt OU 
ae Salt 0 3 
wae 19 3 1D+02 
4.8330D+04 
VW eO0Sop+03 


eo 95 9> 02 
leo D=05 

OS 
-4.8476D-05 

i 4491D-05 
=e 5104 D—- O05 
= aOood p= Os 


THE REDUCED FEEDFORWARD SAT. 


= ee oe On 
120 9290-08 
0.0 
BO 209D—05 
—leo0gcU—05 
1.0794D-04 
26 9207005 


ANALYSIS COMPLETE... YCUR REDUCELE SGisitetei 
IN A FILE NAMED 


NO 


LECLES, 


3.4001D+00 
=2. 20 1o0=01 
050 
= Ole 3 LODO 
I Sls 24 DUS: 
Dios sie Dace 
9.8556 D-04 


(DM) 


1203 open 
3.1109D-—08 
0.0 
-4.4940D-03 
=9253/97 Doce 
-4.0947D-04 
3027 710508 


S23 


1. 7019D+ 00 
4. 8490D-03 
0.0 
-9. 5892p 
=3. I60GD 06 
-7.9878D-04 
-9.4642D-05 


WOPTMA TR SAE awe 


DO YOU DESIRE ANOTHER RUN? 


THE RESULTS OF THIS RUN WILL 


REPLACE YOUR LAST OPTUATR oa 


TYPE eye ci, 


Re T=02 307057732 12 06sec 


RECORD OFF 


END RECORDING OF TERMINAL SESSION 


This 1S a typical computer terminal session in 
the user requested output to 


some of the user protection 


during this session. 


Following this analysis, 
OPTMAT data 


reduced model data. 


renamed 


able characteristics, including stability, 


to the full system eigenvalues. 


raAceoLny 


and further 


and 


the 
OPTS YSxX 


comparison of 
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WNO!M, 


the screen. 


Note also, 


-4.1189D+03 
-1. 3105058 
0.0 
-4.4513D-01 
9.6895D-01 
-3. 1434D-01_ 
-8.3787D- 02 


HAS BEEN SAvae 


wh Zen 
that 


features of OPTKED are involked 


OPTMATR data file was 


was executed 


using 


and 


the 


The eigensystem was analyzed for desir- 
and was compared 
This analysis proved satis- 
the full 


reduced 


models was accomplished by obtaining Bode diagrams for the 
system open loop transfer functions. Figures 4.1 and 4.2 
illustrate the full system and reduced system Bode flots 
respectively. The similarity, especially within the control 
bandwidth, is readily apparent. The frequency responses 
shown are for system input 1 and output 1. Analysis of 
other input/output ccmbinations revealed the same favorable 


comparisons. 
2. Example of a Ninety-eighth Order System 


The X-29A control system includes a ninety-eightn 
order dynamics matrix. The méasurement scaling matrix is 
constructed such that there are two measured outputs. These 
outputs represent the system with and without a notch 
filter, respectively. Full system matrix data as weli as 
eigensysten data was obtained for this model from 
NASA-—Edwards. Full system data was first entered into 
OPTSYSX using subroutine SETUP and the overall system was 
analyzed. Sixty significant states were identified and are 
Shown in the following recorded terminal session. This 
example also illustrates program execution where no "pl" 
matrix iS input and the reduced system states are entered by 
means of file data. During this session the user also 
chooses to have the general reduced model data printed toa 


mrsting file. 


BEGIN RECORDING OF TERMINAL SESSION 
meet —0.01/0.02 12:23:56 
OPTRED 
Poe DEF O5 TERM 
FILEDEF 06 DISK OPTRED LISTING AT 
PEGEDEF O7 DISK STATES DATA Al 
PEG EDEF O09 DISK OPTMAT DATA Al 
FILEDEF 10 DISK OPTMATR DATA A1 
SeeaAL TXATLIB FORTMOD2 MODZEEH IMSLDP NONIMSL 


aD 


LOAD OPTRED (Sia 
EXEGUTLONSE EGE Sr 


OPTRED WILL COMPUTE A REDUCED ORDER MODEL FROM FULL SYSTEM 
UPN wGgH HM) AND "D" (IF INPUT) MATRICES. 


THE FULL SYSTEM MATRICES MUST BE READ FROM AO Fite 
"OPTMAT DATA" AS CREATED BY THE OPTSYSk PROG ae 


YOU MUST ALSO ENTER THE DESIRED REDUCED ORDER (NUSBER OF 33a 
AND THE ACTUAL STATE #'S (IN ASCENDING ORDER) WHICH REPRE 
THE RELUCED MODEL. THE ORDER OF THs aEDUCED JOD i. 

LESS THAN THE ORDER OF THE FULTS Ye 2h 


DO YOU WISH 20 CCNT TNE: 


DYER, VYESYeOn eon 
YES 
DO YOU WISH RESULTS TO SCRED NR Cia. 


NCTE: A DATA FILE CONTAINING THE REDUCED ORDER DATS 
WILL BE GENERATED REGARDLESS OF YOUR ANSWER. 


TYPE "S" FOR SCREEN Of “Dole 


D 
THE ORDERD OF THE SUL SS Pees eee 
THE NUMBER “OF SCONTROLS 20s ae 
THE NUMBER OF OBSERVATIONS 2es: 2 
A "DW MATRIX WiLL NOT BED Pheu 
DO YOU WISH TOeCCNTING 
TYPE OMY ES" SOne NO 
YES 
ENTER THE DESIR&£ED REDUCED CRDER OF THE es 
fs 
60 


DO YOU WISH TO INPUT DESIRED STATES FOR YOUR REDUCED OF DES 
MODEL FROM A DAT Aree: 
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Deere revs eee NAMED "STATES DATA Ai" IN FIXED FORMA 
Pern PORMAT [sea "iso" PER 72 CHARACTER LINE. 
Pee oe OR NOY. 
eS 


Mee kEDUCED MODEL STATES ARES; 

1 Z 3 4 5 6 1z 13 14 15 43 4 
46 47 48 49 50 51 a2 53 54 55 62 63 
65 66 67 68 69 70 71 72 73 74 irs 76 
78 no 80 81 82 83 84 85 86 oF 88 89 
91 92 G3 94 oie) 96 s) aL ks 


WORvOUTWon TO CHANGE ANY OF THE SIGNIFICANT STATES? 


eee pole OR NOt. 
NO 


Preece tolS COMPLETE... YOUR REDUCED SYSTEM DATA HAS BEEN SAVED 
IN A FILE NAMED "OFTMATR DATA". 


DO YOU DESIRE ANCTHER RUN? 


Pao io7e one REOULTS OF THiS RUN WILL 
REEieer YOUR GAS CPINATR DATA. 


Ee eo Ore NOt 
NO 


Mette Oe TREO eto NOW TERM ENAT ED. ..<< 6 « 


eee —17.70/18.51 12: 25:07 
meCORD OFF 
END RECORDING OF TERMINAL SESSICN 


The eigensystems for both systems were analyzed for 
Stability and other desireable characteristics. Bode piots 


were obtained for both input/output combinations and these 
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Te 


45 
64 
vl 
s0) 


graphs for the full and reduced systems .are included as 
Figures 4.3 through 4.10. The input/output combination 
(1/1) shows the response of the system without the notch 
filter while the (1/2) input/output combination depicts the 
freyguency response of the system after the addition of the 
notch filter. These figures clearly indicate excellent 
agreement between the full and reduced order systems in the 


frequency domain. 
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A. CONCLUSIONS 


The computational capabilities of OPTRED were tested by 
analyzing two significant large order systems. The determi- 
nation of the validity of both full order and reduced orcer 
results was greatly simplified by the eigensystem data for 
the X-29A provided by NASA-Edwards. As discussed by Laptas 
Meet. 6 |], OPTSYSX failed to produce valid open-loop transfer 
function zeroes when analyzing very iarge systems. However, 
through modal analysis and with the aid of the NASA data, 
the extraneous zeroes were easily identified and the open 
loop data created by OPTSYSX was corrected. 

The results of OPTSYSX and OPTGRAPH indicate a very 
strong Similarity between the full system and the reduced 
system generated by OPTRED both in stability and in the 
frequency domain. The combined analytic capabilities of 
OPTSYSX, OPTGRAPH and OPTRED previde a powerful tool for the 
controls designer for use in the analysis, development and 


the implementation of complex large order systems. 


B. RECOMMENDATIONS 


Based upon the results of this thesis the following 


areas are recommended for further study and research. 


1. Comparison of System Time Response 


Although the reduced model must closely resemble the 
full system in the frequency domain, it should also compare 
favorably in the time domain. OPTCALC, the time response 
program, requires several modifications to enable its use 


for analysis of systems as large as those presented in this 


4g 


thesis. The addition of this frogram to the analysis tech- 
nique will provide a further means of verification of the 


accuracy and fidelity of the reduced order model. 


2. Interface with 


The Naval Postgraduate Schools Aeronautical 
Engineering Department has developed a collection of anal- 
ysis programs including OPTSYSX, OPTGRAPH and OPTCALC. 
These programs exist ina single "module" which enables the 
user to interactiveiy execute these programs during a single 
computer terminal session. The feasibility of adding OPTRED 
to this module should be investigated. If necessary, these 
programs would require re-dimensioning and could create a 
very large memory size requirement. Additions and modifica- 
tions would also be necessary in the OPTSYSX EXEC to enable 
the execution of OPTRED and then transfer execution to 
OPTSYSX for the analysis of the reduced model. T hate 
controlling EXEC must also manage the data files upon which 


the individual programs depend. 


3. Alternate Order Reduction Technigues 


2 —— ee a oe ce — ee BSS — a co ee a ee SS ee eee a 


Finally, other methods of order reduction should be 
investigated. Comparisons of reduced order models which are 
generated by different techniques would provide the controls 
designer with additional means ci reduced model verification 


and analysis. 
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